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^ ■ Abstract 
(N 

In this paper, we propose a linear precoder for the downlink of a multi-user MIMO system with multiple 
users that potentially act as eavesdroppers. The proposed precoder is based on regularized channel inversion 
(RCI) with a regularization parameter a and power allocation vector chosen in such a way that the achievable 
secrecy sum-rate is maximized. We consider the worst-case scenario for the multi-user MIMO system, where 

>: 

CO . the transmitter assumes users cooperate to eavesdrop on other users. We derive the achievable secrecy sum- 

rate and obtain the closed-form expression for the optimal regularization parameter oils °f tne precoder using 
^ — , ■ large-system analysis. We show that the RCI precoder with oils outperforms several other linear precoding 

o: 

, schemes, and it achieves a secrecy sum-rate that has same scaling factor as the sum-rate achieved by the 

optimum RCI precoder without secrecy requirements. We propose a power allocation algorithm to maximize 
the secrecy sum-rate for fixed a. We then extend our algorithm to maximize the secrecy sum-rate by jointly 
optimizing a and the power allocation vector. The jointly optimized precoder outperforms RCI with «ls an d 
equal power allocation by up to 20 percent at practical values of the signal-to-noise ratio and for 4 users and 
4 transmit antennas. 
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I. Introduction 

In current practical multi-user MIMO systems such as LTE and 802.1 In, securing transmitted data 
from nearby eavesdroppers is critical. In these systems, security is achieved using potentially vulnerable 
network layer cryptography techniques. The vulnerability is due to a reliance on the limited resources of 
the eavesdropper and on the unproven computational complexity of inverting the encryption algorithms 
flU. To enhance the protection of transmitted data and achieve perfect secrecy, methods exploiting the 
channel, known as physical layer security, have been proposed. 

Physical layer security techniques were proposed to protect the data from eavesdroppers for several 
network topologies in [|2l- [fl4ll . In [O, 0, a three-terminal network consisting of a transmitter, an 
intended user and an eavesdropper, known as the wiretap channel, was considered. The authors derived 
the secrecy capacity, where the message is transmitted reliably to the intended user while the rate of 
information leakage to the eavesdropper vanishes asymptotically with the code length. The secrecy 
capacity of the wiretap channel was derived for the MIMO case in JH-|l6l when all terminals had 
full channel state information. It was shown in Q-QT]] that the transmission of artificial noise, as 
well as adaptive encoding, is an effective method to reduce the eavesdropper's signal-to-noise ratio 
when the eavesdropper's channel is not known by the transmitter. Recently, physical layer security was 
also extended to multiuser networks where the eavesdropper is not an intended user [[121 . [fT3l and to 
two-user networks where the intended users are also eavesdroppers [|T4l . The secrecy capacity region 
for multi-user networks where any number of intended users are potentially eavesdropping remains 
an open problem. Moreover, the achievable secrecy rates of such multi-user networks with practical 
transmission schemes are also unknown. 

Suboptimal precoding schemes have proven to be practical and effective in controlling interuser 
interference for the downlink of multi-user MIMO networks [fT51 - [|23l . While the sum-capacity of 
multi-user MIMO networks without eavesdroppers is achieved using dirty paper coding E4l . it requires 
high-complexity coding schemes [|25l . Linear precoding schemes were proposed as a low-complexity 
alternative for multi-user MIMO downlink implementations [26]. A popular and practical linear pre- 
coding scheme to control interuser interference is channel inversion (CI) precoding, sometimes known 
as zero forcing precoding [fT5l , [fTBI . To increase the sum-rate performance of the CI precoder, the 
regularized channel inversion (RCI) precoder was proposed to tradeoff the interuser interference and the 
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desired signal through a regularization parameter [fTTll . Linear precoding schemes were also proposed 
to achieve secrecy in single-user MIMO networks Hl-Q. 

In ll27l . the use of linear precoding was proposed to achieve physical layer security in a multi-user 
MEMO system. For RCI precoding, the authors obtained an achievable secrecy sum-rate as a function 
of the singular values of the channel. A bound on the optimal regularization parameter of the precoder 
was also given in the large-system regime. 

In this paper, we consider the multi-user MIMO downlink with multiple single- antenna users that 
cooperate and jointly eavesdrop on other users, and we propose a linear precoder based on RCI. We 
use large-system analysis with an approach different from that of ||27ll , and we derive the optimal 
regularization parameter ccls and the corresponding achievable secrecy sum-rate. Numerical results 
confirm the accuracy of the large-system analysis, even when applied to a number of users as low as 
4. Moreover, the RCI precoder with «ls outperforms several other linear precoding schemes. In fact, 
it achieves a secrecy sum-rate that has same scaling factor as the sum-rate achieved by the optimum 
RCI precoder without secrecy requirements. We then propose an iterative power allocation algorithm 
to obtain the maximum secrecy sum-rate for fixed a. We extend our algorithm to maximize the secrecy 
sum-rate by jointly optimizing the regularization parameter a and the power allocation vector. The 
proposed power allocation algorithm outperforms RCI with a LS and equal power allocation (RCI-EP) 
by up to 20 percent at practical values of the SNR and for 4 users and 4 transmit antennas. 

Throughout the paper we use the following notation: bold uppercase (lowercase) letters denote 
matrices (column vectors); (-) T and denote matrix transpose and conjugate transpose, respectively; 
the trace of a matrix is denoted by tr{-}, and the Euclidean norm of a vector is indicated by || ■ ||; 
E[-] denotes the expected value of the random variable in the brackets, CJ\f(fi, & 2 ) denotes circularly 
symmetric complex-Gaussian distribution with mean // and variance a 2 , and we use the notation 
[■]+ = max(-,0). 

II. System Model 

We consider the downlink of a narrowband multi-user MIMO system, consisting of a base station 
(BS) with M antennas which simultaneously transmits K independent confidential messages to K 
spatially dispersed single-antenna users. Transmission takes place over a block fading channel, where 
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the coherence time of the channel is much longer than one symbol interval. In this model, the 
transmitted signal is x = [xi, . . . , xm] T G <C Mx1 , and the received signal at user k is given by 

M 

Vk = y^ y h ktj Xj + n k (1) 

where h k j ~ £/V(0, 1) is the fading gain between the j-th transmit antenna element and the A;-th user, 
and n k ~ £/V(0, a 2 ) is the noise seen at the /c-th receiver. The corresponding vector equation is 

y = Hx + n (2) 

where H = [h k j] is the K x M channel matrix, y = [yi, . . . , y K ] T and n = [ni, . . . , n^-] T . We impose 
the long term power constraint E[||x|| 2 ] = 1, assume that E [nnt] = <j 2 l, and define the SNR p = 1/a 2 . 
The transmitted signal x is obtained at the BS by performing a linear processing on the confidential 
messages u k , k — 1, . . . , K. 

It is required that the BS securely transmits each confidential message u k , ensuring that the unin- 
tended users receive no information. This is performed at the secrecy rate R s ,k, defined as follows. Let 
Pr(£ n ) be the probability of error at the intended user, m be a confidential message, y™ be the vector 
of all signals received by the eavesdroppers, and if(m|y") be the corresponding equivocation. Then 
a (weak) secrecy rate R 3 k for the intended user is achievable if there exists a sequence of (2 nHs fc , n) 
codes such that Pr(£ n ) — > and ^H(m\y™) < ^H(m) — e n with e n approaching zero as n — > oo O. 

In general, the behavior of the users cannot be determined by the BS. As a worst-case scenario, in our 
system we assume that for each intended receiver k the remaining K — 1 users can cooperate to jointly 
eavesdrop on the message u k . For each user k, the alliance of the K — 1 cooperating eavesdroppers 
is equivalent to a single eavesdropper with K — 1 receive antennas, which is denoted by k. 

III. Linear Precoding 

In this section, we derive an achievable secrecy sum-rate for the multi-user MIMO downlink with 
malicious users by using a linear precoder. Although suboptimal, linear precoding schemes are of 
particular interest because of their low-complexity implementations and because they can control the 
amount of crosstalk between the users [fT51 - lfT8l . We then specialize and obtain the secrecy sum-rate 
achievable by the RCI precoder. RCI is a linear precoding scheme that was proposed to serve multiple 
users in the multiuser MIMO downlink channel, which has better performance than plain channel 
inversion, especially at low SNR IfTTl . 



A. Preliminaries 

In linear precoding, the transmitted vector x is derived from the vector containing the confidential 
messages u = [u%, . . . , uk] T through a deterministic linear transformation (precoding) |[T5l - [fT8l . We 
assume that the entries of u are chosen independently, satisfying E[|ufc| 2 ] = 1. We assume spatially 
homogeneous users, i.e. each user experiences the same received signal power on average, thus the 
model assumes that their distances from the transmitter are similar. 

Let W = [wi, . . . , wjf] be the M x K precoding matrix, where is the A;-th column of W. Then 

the transmitted signal and the power constraint are, respectively: 

1 1 K 

x = — Wu = — — ) w fe w fc , (3) 

E[||x|| 2 ]=lE[||Wu|| 2 ]=lf;||w fc || 2 = l, (4) 

^ ^ k=i 

where 7 = tr{W^W} is the long-term power normalization constant. 

B. Achievable Secrecy Sum-Rates with Linear Precoding 

By employing the linear precoding in ©, the signals observed at receivers k and k are, respectively 

Vk = -^=h[w k u k + — — h[w jUj + n k 

^ (5) 



-Eh, 



where = [ni, . . . , n k -i, Uk+i, ■ ■ ■ , uk] , h^. is the k-th row of H, and is a matrix obtained from 
H by eliminating the A;-th row. The channel in © is a multi-input, single-output, multi-eavesdropper 
(MISOME) wiretap channel [6J. The transmitter, the intended receiver and the eavesdropper of this 
MISOME wiretap channel are equipped with M, 1 and K — 1 virtual antennas, respectively. Due to 
the simultaneous transmission of the K messages, user k experiences noise and interference from all 
the Uj, j ^ k. 

In the following, we derive an achievable secrecy sum-rate R s for the multi-user MIMO system 
with malicious users. Although the design of codes for the multiuser MIMO channel with security 
constraints is not the focus of this paper, we prove the achievability of R s with a code construction 
based on independent codebooks and linear precoding. 



Lemma 1 (Codebook construction): An achievable secrecy sum-rate R s for the multi-user MIMO 
system with malicious users is given by 



K 

A 



R s — R s ,k, (6) 



k=l 



where R Sjk is an achievable secrecy rate for the A;-th MISOME wiretap channel ©, k = 1, . . . , K. 

Proof: Assume that the BS uses independent codebooks for each user, where each codebook is 
a code for the scalar wiretap channel J6]|. The confidential message u k is obtained as a codeword 
independently drawn from the code C k , corresponding to the A;-th user. The rate R s ^ of the code C k is 
chosen according to the secrecy rate achievable for user k in the presence of eavesdropper k, i.e. by the 
secrecy rate achievable for the MISOME wiretap channel ©. The existence of such code is guaranteed 
by the definition of secrecy rate 01 . To construct the vector codeword for the broadcast channel, the 
scalar codewords for each MISOME wiretap channel are stacked according to u = [m, . . . , %] r , and 
no additional binning is required. The vector u is then linearly precoded as in ©, which means that 
each message u k is transmitted by beamforming (i.e. signaling with rank one covariance) along the 
direction w k . The secrecy sum-rate R s is then by definition the sum of the simultaneously achievable 
secrecy rates R S:k . ■ 
Lemma 2: An achievable secrecy rate R s * for the MISOME wiretap channel © is given by 



R 



s . k 



log 2 (l + SINR fe ) - log 2 (l + SINR^ 



(7) 



where SINR^ and SINR^ are the signal-to-interference-plus-noise ratios for the message Uk at the 
intended receiver k and the eavesdropper k, respectively. 

Proof: By noting that the MISOME wiretap channel © is a nondegraded broadcast channel 0, 
the secrecy capacity is given by J3J: 

C s = max I[u k ]y k ) - l(u k ]yr) (8) 

where I(x; y) denotes mutual information between two random variables x and y. The secrecy capacity 
C s is given by the difference of the mutual informations at the intended user and at the eavesdropper, 
respectively. C s is achieved by maximizing over all joint probability distributions such that a Markov 
chain u k — > w k u k — » yk,y k is formed, where u k is an auxiliary input variable. By evaluating ([8]) with 



u k ~ CAf(0, 1) and with the linearly precoded data w k u k , we obtain an achievable secrecy rate R Sjk 
for the MISOME wiretap channel © given by 



R s ,k — 



(9) 



where (a) follows from w k u k being a deterministic function of u k (6l. Equation © then follows from 
© and from the statistics of u k . ■ 

From equation © it is clearly observed that for high-performance linear precoder design an efficient 
tradeoff between maximizing SINR fc and minimizing SINR^ is required. 

Theorem 1: A secrecy sum-rate achievable by multi-user MIMO linear precoding is given by 



A" 



k=l 



log 2 1 + 



h fc w fc 



7^ 2 + Ej^fc 



- log 2 1 



| H fe w fc| 

7<T 2 



(10) 



Proof: By using Lemma [Hand Lemma [2l we have that an achievable secrecy sum-rate is obtained 
as the sum of the secrecy rates R Sjk in ©. A lower bound on the quantities R s k can be obtained by 
considering a genie-aided eavesdropper which observes not only the signals received by its K — 1 
antennas, but also all the confidential messages Uj,j ^ k. Such channel clearly has an achievable 
secrecy rate smaller than the original channel. The genie-aided eavesdropper k can perform interference 
cancellation, and it does not see any undesired signal term apart from the received noise n^. 

According to the previous considerations, the signals at the intended receiver and the eavesdropper 
of the k-th equivalent MISOME wiretap channel become, respectively: 

— h[w k u k + — = 2J hlw jUj + 



Vk 



7 



7 j& 



n k 



(11) 



y~k = -^H^w fc u fe + 



For the fc-th equivalent MISOME wiretap channel in (fTTj) . the SINRs at the intended user and the 
eavesdropper are, respectively: 

t 2 



SINR A 



7^ 2 + E, ¥fe 



2 • 



SINRfc 



|H^w fc | 
7<j 2 



(12) 



(13) 



Since the noise in in (fTTI) is spatially white, the optimal receive filter at k is the matched filter 
(H^w fc )^. Equation (fT3T ) then follows. For a given channel H, substituting (fT2l) and (fT3T ) into (Q and 
then into © yields ([TO)) . ■ 
For the remainder of the paper, we refer to equation (fTOl) as the secrecy sum-rate. We note that it 
depends on the choice of the precoding matrix W, as well as on the channel H and the noise variance 
a 2 . A possible choice for W, based on regularized channel inversion, is discussed in the following. 

C. Achievable Secrecy Sum-Rates with Regularized Channel Inversion 

We now consider RCI precoding for the multi-user MIMO downlink with malicious users. Although 
CI precoding can achieve secrecy by canceling all signals leaked at the unintended users, this comes 
at the cost of a poor sum-rate. The RCI precoder has better performance than plain CI, particularly 
at low SNR IfTTl . For each message RCI precoding achieves a tradeoff between the signal power 
at the A;-th intended user and the crosstalk at the other (K — 1) unintended users for each signal. The 
crosstalk causes interference to the unintended users. In the case when the unintended users are acting 
maliciously, the crosstalk also causes information leakage. Therefore, RCI achieves a tradeoff between 
signal power, interference, and information leakage. 

With RCI precoding, linear processing exploiting regularization is applied to the vector of messages 
u IfTTl . The RCI precoding matrix is given by 

W = H f (HH f + aI K ) -1 . (14) 

The transmitted signal x after RCI precoding can be written as 



x = — Wu = — H f (HH f + aI K ) 1 u = — (H f H + aI K ) 1 H f u. 
The latter passes through the channel, producing the vector of received signals 



(15) 



y = — H (H f H + aI K ) 1 H f u + n. (16) 
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The function of the real nonnegative regularization parameter a is to improve the behavior of the 
inverse, although it also produces non-zero crosstalk terms in (fT6l) . 
Using RCI precoding, the SINRs (fT2l) and (TT~3T) at the intended user k and the eavesdropper k become 

h£ (HtH + «I^) _1 h A 

SINR fc 



7^ 2 + Ej^fc 



2> (17) 
hi (HtH + aK-r 1 h, 



SINR, 

where 



H^HtH + aljr) 1 h k 



(18) 



7 = tr{H t H(H t H + o;I^) 2 } . (19) 

To simplify (fT71) and (fl"8~l) . we introduce the quantities 

A k = h{ (nlU~ k + al K ^j ~* h fe and (20) 

Bk = ^ (h+H s + aI K ) ~ l H~H^ (h£h s + ol*) ' h*. (21) 

It is then possible to express (fTTT) as [|28l 

SINR fc = h - (22) 

Bk + ya* (1 + A k f 



In a similar fashion, we rewrite (1181) as 

SINRr = — J. (23) 

By substituting (l22~l) and (1231 into © and then into © we obtain the following expression for the 
secrecy sum-rates achievable with RCI precoding 



K r 



fc=i 



logo 1 + ~ o - logo 1 -r o 



IV. Large-System Analysis 



(24) 



In this section, we consider the performance of the RCI precoder in the large-system regime, where 
both the number of transmit antennas M and the number of receivers K approach infinity in a fixed 
ratio. We derive closed-form expressions for the optimal regularization parameter and the optimal 
secrecy sum-rate achievable with RCI in the large-system regime. We then compare the secrecy sum- 
rate achieved by the optimized RCI precoder to several other linear precoding schemes. Finally, we 
evaluate the sum-rate loss due to the secrecy requirements. In this section, we focus on the case 
K = M because the analysis for this case is more tractable, and it is considered to be most important 

ma. 



A. Secrecy Sum-Rates in the Large-System Regime 

We define £ = a/K as the normalized regularization parameter, and note that as K 
quantities (fT9l) . (l20l) and (ED converge (almost surely) to J28l 

hm 7 = ^(0+^(0, 



lim A fc = # (f ) 

if— »0O 



d 



hm 5, = 9 (0 + Z-7Z9 (0 



respectively, where 



<K0 



1 --I 

2 V + £ 2' 



oo, the 

(25) 
(26) 
(27) 

(28) 



By substituting the above expressions in (|22|) and (1231) . one can conclude that as K — > oo, all the SINRs 
at the intended user k and at the eavesdropper k converge to a non-random function of the parameter 
£ and the noise variance a 2 . Moreover, these quantities are the same for all confidential messages Uk, 
as K — > oo. Hence, in the large-system regime, it is possible to write the secrecy sum-rate with RCI 
precoder as 



R, ^ Rs 



A 



K 



1 + 



log 2 



P9(0 2 



1 + 



(i+s(or 



as K — > oo. 



(29) 



B. Selection of the Optimal Regularization Parameter 

The value of the asymptotic secrecy sum-rate R Sy00 in (l29l) depends on the normalized regularization 
parameter £. We now derive the optimal value £ opt that maximizes R s ,oo- 
Lemma 3: The optimal normalized regularization parameter in the large-system regime is given by 

1 



6 



opt 



3p + 1 + v^+T 



(30) 



Proof: The value of £ opt is obtained as the stationary point of the secrecy sum-rate R Sy00 , which 
can be found by setting to zero the derivative of the logarithm in d29l) . by applying some algebraic 
manipulations, and showing that the maximum value is nonnegative. ■ 
As in the case with no secrecy requirements, the value of £ opt is a function of the SNR p. In a 
multi-user channel without secrecy requirements, the choice £ = 1/p is optimal for large K, as it 
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maximizes the sum-rate of the system ifTTll . However, this value is no longer optimal in a multi-user 
channel with malicious users. In fact, because of the secrecy requirements the crosstalk terms appear 
twice in the secrecy sum-rate expression (flOl) . As a consequence, £ = 1/p is too large and gives too 
much crosstalk to the other users. This was proven in [|27l and it is also easily confirmed by the 
following inequality: 

Wt<^<- Vp. (31) 
3p p 

Similarly to the case with no secrecy requirements, £ opt decreases as we increase p. The high-SNR 
asymptote of £ opt is given by 

£ op t ^ ■=-, as p ->■ oo (32) 
3p 

and £ opt tends to zero if p — > oo. 

Unlike the case with no confidentiality, the optimum normalized regularization parameter is upper 
bounded and it does not tend to infinity as p tends to zero. The low-SNR asymptote of £ opt is 

£ op t = ~, for p = 0. (33) 

In the remainder of the paper, we will denote by ccls = K£ opt the unnormalized large-system 
regularization parameter. 

C. Optimal Secrecy Sum-Rate 

It is now possible to obtain an expression for the optimal secrecy sum-rate of the RCI precoder in 
the large-system regime. The optimal secrecy sum-rate is a function of the SNR p and the number of 
users K only. 

Theorem 2: The optimal secrecy sum-rate Rf^ achievable by the RCI precoder in the large-system 
regime is given by 



K s,oo - maxK Si00 - A log 2 4 (4p + 1) ' 



Proof: Equation (1341) is obtained by substituting (1301) in (1291 ) and applying some algebraic ma- 
nipulations. ■ 
The secrecy sum-rate -R^ 1 in (|34l) satisfies 

Rf™>0 Vp>0, (35) 
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and the high-SNR asymptote of -R^ 1 is given by 

K 27 K 

R fS ~ y lo S2 ^ + y lo S2 P, as p ^ oo. (36) 

Therefore, in the large-system regime the secrecy sum-rate for optimal £ scales logarithmically with 
high SNR, and it scales linearly as K/2 with the number of users. 

We have shown that although for K — > oo the number of eavesdroppers K—l for each message tends 
to infinity, a positive secrecy sum-rate is still achievable. This occurs because the number of transmit 
antennas M = K also tends to infinity, and it is larger than the number of eavesdroppers. Therefore, 
for each message the transmitter is able to control the amount of interference and information leakage. 

We now compare the per-user secrecy rate achieved by RCI to the secrecy capacity of the MISOME 
channel, C^isome^ m ^ high-SNR regime. The former is obtained by dividing (|36l) by the number 
of users K, and it can be further approximated by 

dRCI 

—j-r- ^ ^ l0g2 aS P ~^ 00 • ^ 37 ^ 

The value of (7MISOME was obtained in J6]|, and it can be approximated by the following lower bound 

Cf ISOME >~log 2 P, asp ^00, (38) 

which is tight at high SNR (6]|. We note that in C*misome f rom [jgj a single-user system is considered. 
Therefore, only one message is transmitted to one legitimate user, and the user does not experience 
any interference. For large SNR, the RCI precoder achieves a per-user secrecy rate which is the same 
as the secrecy capacity of a single-user system. 

D. Comparison to Other Linear Schemes 

In the following, we compare the secrecy sum-rate in (1341) achieved by the RCI precoder with £ opt to 
the secrecy sum-rates obtained from (|2"91) by using: 1) £ = (CI precoder), 2) £ — > 00 (matched- filter 
precoder) and 3) £ = 1/p (optimum RCI precoder without secrecy requirements). 

The aim of the CI precoder is to cancel all the interference and information leakage, therefore 
yielding to a secrecy sum-rate that coincides with the sum-rate. We note that for the CI precoder it is 
£ = 0, and the precoding matrix is given by 

W = H t (HH t ) _1 . (39) 
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In order for the inverse in d39l) to exist, it is required that K < M. 

The secrecy sum-rate achieved by CI in the large-system regime grows at most sublinearly with 
K — > oo. In fact, 

lim lim — = 0. (40) 

f-K)iC-Kx> K 

This result is consistent with ll28l . where it was shown that the CI precoder performs poorly in the 
large-system regime when the number of antennas equals the number of users. 

Similarly, we calculate the secrecy sum-rate achieved when £ — > oo (matched- filter precoding). Here, 
the transmitter beamforms in a direction such as to maximize the signal strenght of each user, without 
taking into account the interference it creates and the amount of resulting information leakage. The 

secrecy sum-rate achieved by matched-filter precoding in the large-system regime is zero. In fact, 

2p+ 1 



lim lim — y 



logo x 

. (P+l) 2 



0. (41) 



Clearly, matched-filter precoding performs poorly compared to the optimal RCI precoder. This is due 
to the intended user suffering from a large amount of interference, while the eavesdroppers may cancel 
the interference by cooperating. 

Finally, we consider £ = which is the value that maximizes the sum-rate of the system without 
secrecy requirements IfTTl . The secrecy sum-rate R° s oo achieved by RCI with £ = 1/p in the large- 
system regime is given by 



«..oo = K\0 g2 ^Tl) • (42) 

We observe that the RCI scheme with £ = 1/p outperforms the CI and the matched-filter precoding 
schemes in the large-system regime, but it is suboptimal compared to the use of £ opt . For high SNR, 
the per-antenna secrecy sum-rate gain provided by using £ = £ opt in place of £ = 1/p is given by 

lim R ^- R °^ = log2 b/1 ~ o.38 bits. (43) 

p->oo K 4 

E. Secrecy Loss 

We now consider the secrecy loss, i.e. the sum-rate loss due to the secrecy requirements. We define 
this as the difference between the optimal sum-rate R ^ without secrecy requirements and the secrecy 
sum-rate -R^ 1 in (|34l . The sum-rate is obtained with RCI and £ = 1/p, and it is given by ll28l 



^ = /rlog2 i±V^±l. (44) 
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The high-SNR asymptote of the sum-rate in (|44|) is 



2 

For high SNR, the per-antenna secrecy loss is given by 



Kc - 7T log 2 P, as p ^ oo. (45) 



^rci \ g4 

lim °° T ^ = - log 2 — « 0.62 bits. (46) 

p^oo K 2 62 27 

By comparing (1431) to (|36l ), one can conclude that the secrecy requirements do not change the linear 
scaling factor for large SNR. In other words, the RCI precoder with £ opt achieves a secrecy sum-rate 
that has same scaling factor K/2 as the sum-rate achieved by the optimum RCI precoder without 
secrecy requirements in ifTTI . The RCI precoder with £ opt can achieve secrecy with a penalty in terms 
of the per-antenna sum-rate given by (|46l ). The secrecy loss (|46l corresponds to a power loss of a 
factor 64/27 ~ 3.75dB. Therefore, the RCI precoder with £ opt can achieve secrecy without reducing 
the sum-rate of the system, as long as the transmitted power is increased by 3.75dB. 

V. Power Allocation 

In this section, we consider power allocation for the RCI precoder. We first propose a new algorithm 
to obtain the power allocation vector p which achieves the optimal secrecy sum-rate with a fixed 
regularization parameter a. We then extend our algorithm to jointly optimize p and a. 

A. Achievable Secrecy Sum-Rates 
We consider the RCI precoding matrix with arbitrary power allocation given by 

W p = WD = H^HH* + aI) -1 D, (47) 

where D = diag( v /p), and p = [pi, . . . ,Pk] T is the power allocation vector. The vector p must be 
chosen such that the power constraint tr{W p W p } = 1 is met. Clearly, (|47T) generalizes the RCI 
precoder W with equal power allocation (RCI-EP) in (TT4l) . 
When the precoder W p is used, the SINR at the k-th intended user, given by (fl"2l) . becomes 



SINR, = - Pfc| ^ Wfc|2 (48) 



and the SINR at the eavesdropper k, given by (fl"3l) . becomes 



SINR E = »'' H f'»' = P'S*»l h WP (49) 



a" a* 
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From (1481) and (1491) , we obtain the achievable secrecy sum-rate with power allocation 

K 



k=l 



(50) 



5. Power Control 

To obtain the optimal power allocation vector p, we are required to solve the non-convex optimization 
problem 

maximize R^ A (p) 

(51) 

subject to tr {W p W p } < 1, 
where i?^ A (p) is given by (1501) , W p is given by (1471) , and the maximum total transmit power over 
all antennas is one. In the following, we will ignore the notation [■] + in (l50l) in the maximization 
problem. In fact, any negative term in the sum can be replaced by zero (thus increasing the sum) by 
using p k = which is always feasible. 

We now reformulate the problem (|5TT) by applying the transformation p k = \ogp k , k — 1 . . . , K, 
and obtain the optimization problem 

maximize R^ A (p) 



(52) 



subject to tr {W p W p } < 1, 

where p= [pi,... ,Pk] t - 
Lemma 4: The second term of the objective function, i?^ A (p), of (l52l) is concave. 
Proof: The second term and its first and second derivatives are 



-log, (1+SIN%) = -log 2 (^l+ ePfcE ^ 2 |hjWfc 
dlog 2 (1 + SIN%) _ e& |hjw fc | 2 



d 2 log 2 (1 + SIN%) _ e& £^ fc V 



lo S2 e, (53) 



- / x2 -log 2 e<0. 



9pi (^ 2 + e ?fe E j¥ Jh]w fe | 2 

Hence, by the second order condition [30, §3.4.3] , — log 2 (l + SINR^) is concave 



In order to solve the problem (1521) , we consider a modified version of the method as in [29] and 
[1301 which is based on a reformulation of (1521 . This approach guarantees an improvement in the 
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performance over the standard high-SNR approximation in fading channels 11301 . In order to obtain 
the reformulation, we use the following bound obtained in If29l 

alogz + b< log(l + z), a= Z ° and b = log(l + z Q ) - - — logz , (54) 

1 + z 1 + z 

for some z > 0, with equality when z = z . 

Lemma 5: With the change of variables pk = \ogpk, k = 1 . . . , K, the lower bound 

Gfc io g L e ": |hlwfcl ,! , VA<io S Ji+ ^ Pfc|h K j , (55) 



(56) 



lQ g2 VE^ fc ^|hiw,P + o*J log2 - ^ E^Piihlwil 2 + o 3 
is concave in pk, k = 1, . . . , K. 

Proof: The result follows immediately using the method in Lemma @] ■ 
We showed in Lemma @] that the second term of (l50l) is concave by the second order condition. By 
using the lower bound in (1551 ) for the first term of (1501 ), we obtain a concave objective function. Since 
the constraints are affine, the optimization problem arising from (|52|) and the bound (1551) is a convex 
optimization problem. This convex optimization problem is given by 

.. ^fa*. ( e^|hjw fc | 2 \ 6 fc / ^E^Jhjwfc 
maximize > log — +- log 9 lH J ^ - — 

subject to tr {Wj,W p } < 1 
The power allocation vector can then be obtained using Algorithm 1 in Table HI To show that 
Algorithm 1 converges monotonically to a local optimum, we note the constraint is the same for both 
the i-th and (£ + l)-th subproblems. Hence, the solution of the t-th subproblem (|56l) is also feasible for 
the (t+l)-th subproblem (l56l) . Moreover, by the bound in (l54l) . the objective function is monotonically 
increasing and converges to a local optimum. 

C. Proposed Precoding Scheme 

Having established an algorithm to determine the optimal power allocation vector p for a fixed a, 
we now obtain our precoding scheme by considering the joint optimization of a and p. The joint 
optimization problem can be written as 

maximize i?f A (p, a) 

P,a ' (57) 

subject to tr {WjW p } < 1. 
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Even after using the transformation pj. = logpfc, k = 1, . . . , K, the problem (|57|) is non-convex. To 
solve this problem, we propose Algorithm 2 in Table [0 

At each iteration, Algorithm 2 optimizes the regularization parameter a and subsequently the power 
allocation vector p. It is straightforward to prove that Algorithm 2 converges monotonically and it thus 
provides with a locally optimal pair (a, p) for the proposed linear precoder. In Section VI we show 
via simulations that the proposed precoder with jointly optimal regularization parameter and power 
allocation vector outperforms RCI precoding with a L s and equal power allocation (RCI-EP). 

VI. Numerical Results 

In this section, we show the performance of our proposed precoding scheme via simulations. We also 
consider the finite user scenario to show that many results from the large-system analysis in Section 
IV hold for a small number of users. The precoding matrix W was normalized by ^/j, as in (fl9l) . in 
order to meet the power constraint in (HJ). This corresponds to a long-term power constraint, which 
does not require the receivers to know the instantaneous value of 7 [17]. In the following, we denote 
by «ls — -^£opt the large-system regularization parameter, obtained from (|3Q|) . 

Fig. CD compares the large-system regularization parameter a^s to the optimal regularization parameter 
aps f° r a finite number of users. The value of aps was found by using single- variable numerical 
optimization to maximize the mean value of the secrecy sum-rate in (1241) . The figure shows the finite- 
system and large-system regularization parameters at practical SNR values for four different numbers 
of users: 4, 8, 16 and 32. We observe that as the number of users K increases, the value of «fs 
approaches the large-system regularization parameter axs- 

In Fig. |2] we demonstrate that using the large-system regularization parameter «ls in a finite-size 
system does not cause a significant loss in the secrecy sum-rate compared to using a regularization 
parameter ctFs(H) optimized for each channel realization. The figure shows the complementary cu- 
mulative distribution function (CCDF) of the normalized secrecy sum-rate difference between using 
«ls and ccfs(H) as the regularization parameter of the RCI precoder for K = 4, 8, 16, 32 users at 
an SNR of lOdB. The difference is normalized by dividing by the secrecy sum-rate of the precoder 
that uses otps(H). We observe that the average normalized secrecy sum-rate difference is less than 2.4 
percent for all values of K, As a result, the large-system regularization parameter ols may be used 
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instead of the finite-system regularization parameter with only a small loss of performance. Moreover, 
the value of a LS does not need to be calculated for each channel realization. 

Fig. [3] compares the analytical secrecy sum-rate of the RCI precoder in (|34|) to the simulated secrecy 
sum-rate of the RCI precoder with a finite number of users, which is averaged over 10 3 channels. The 
RCI precoder with a finite number of users was obtained by using the regularization parameter c*fs, 
found by simulation, that maximizes the average secrecy sum-rate. We observe that the large-system 
analysis is accurate at low SNR for all values of K. Moreover as K increases, the large-system analysis 
is accurate for larger values of the SNR. 

In Fig. |4] we compare the simulated secrecy sum-rate of the RCI precoder using the large-system 
regularization parameter o;ls with CI precoding lfi~5l and RCI precoding with a = K/p, which 
maximizes the sum-rate without secrecy IfTTl . The sum-rate of the optimal RCI precoder without 
secrecy requirements is also plotted. The figure shows plots for K = 4, 8, 16, 32. We observe that 
CI precoding exhibits a large performance loss compared to the secrecy sum-rate of the optimal RCI 
precoder for large values of K. The RCI precoder with a = K/p outperforms CI precoding, but it is 
suboptimal compared to the RCI precoder that uses ols- We note that although CI precoding achieves 
secrecy in a simple way by completely canceling the information leakage, this comes at the cost of a 
poor sum-rate. Secrecy can be achieved with a significantly larger sum-rate by using the RCI precoder 
with axs. We also observe that the secrecy loss between the sum-rate of the RCI precoder without 
secrecy and the secrecy sum-rate of the RCI precoder is almost constant at high SNR for large K. 
This confirms the result we derived in ( |46l ). Moreover, the value of the simulated per-antenna secrecy 
loss is 0.59 bits for K = 32 and p = 25dB; close to the 0.62 bits suggested by the analysis in (|46|) . 

In Fig. [5] we compare the simulated per-user secrecy rate of the proposed precoder with jointly 
optimized regularization parameter a and power allocation vector p to the RCI precoder with «ls 
and p pt, and to the RCI-EP precoder with a L s- We observe that there is a negligible performance 
difference between the proposed precoder and the RCI precoder with a^s an d p op t- As a result, a low- 
complexity, near-optimal RCI precoder may be implemented by using a LS and optimizing the power 
vector separately. The figure shows that for K = 4, the proposed power allocation scheme always 
outperforms the RCI-EP precoder with «ls by up to 20 percent, and the gain does not vanish at high 
SNR. This occurs because at high SNR £ opt — > and the RCI precoder behaves as a CI precoder, 
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for which the optimal power allocation is waterfilling lfT5l . Hence, equal power allocation for RCI is 
suboptimal at high SNR. Fig. [5] also shows that the proposed power allocation scheme reduces the 
sum-rate loss due to the secrecy requirements. For p > 15dB, RCI with power allocation achieves 
a per-user secrecy rate which is even higher than the per-user rate achieved by the optimal RCI-EP 
without secrecy requirements. Furthermore, Fig. [5] shows the simulated secrecy capacity (7^isome Q f 
a MISOME channel with the same per-message transmitted power. Although (7^isome - ls obtained 
in a single-user and interference-free system [6], at high SNR, RCI with power allocation achieves a 
per-user secrecy rate as large as cmisome^ 

VII. Conclusions 

In this paper, we considered the problem of secret communication in a multi-user MIMO system 
with malicious users. We proposed a linear precoder based on regularized channel inversion (RCI) 
with a regularization parameter and power allocation vector that maximize the achievable secrecy 
sum-rate. The analysis presented in the paper, as well as the simulation results, showed that RCI with 
equal power allocation (RCI-EP) and with the optimal regularization parameter outperforms several 
other linear precoding schemes. Moreover, it achieves a sum-rate that has same scaling factor as the 
sum-rate of the optimum RCI precoder without secrecy requirements. The secrecy requirements result 
in a loss in terms of the sum-rate. This loss can be compensated by the proposed power allocation 
scheme, which increases the secrecy sum-rate compared to RCI-ER 

Part of the analysis presented in this paper focused on the case when the number of users K equals 
the number of transmit antennas M. Generalizing the results to the case when K and M can take any 
value is part of our ongoing research effort. When K > M, the secrecy sum-rate degrades due to the 
increased interference and information leakage. Therefore, it can be useful to design a user scheduling 
algorithm that selects a subset of the users to be served, thus increasing their SINR. However, user 
scheduling cannot reduce the number of potentially malicious receivers in the network, since discarded 
users are still able to eavesdrop. The transmission of artificial noise can limit the eavesdropping ability 
of the discarded users, but it must be harmless to intended receivers. 

Throughout the paper, we focused on the worst-case scenario when the transmitter assumes that 
users cooperate and jointly eavesdrop on other users. We are interested in this scenario because the 



19 



transmitter is unable to predict whether the users are eavesdropping or not. Possible extensions of this 
work include considering a scenario where only some of the users eavesdrop on other users, or where 
users can individually eavesdrop, but without cooperating. In this case, the secrecy rate for each user 
is limited by the eavesdropper that receives the largest information leakage. We leave the analysis of 
these aspects for future work. 
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Fig. 1. Comparison between the large-system regularization parameter q?ls and the value qfs that maximizes the average secrecy 
sum-rate for finite K. 
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Fig. 2. Complementary cumulative distribution function (CCDF) of the normalized secrecy sum-rate difference between using a-ps (H) 
and qls with p = lOdB. 
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Fig. 3. Comparison between the secrecy sum-rate in the large-system regime d34t and the simulated secrecy sum-rate for finite K. 
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Fig. 4. Comparison between the RCI precoder with qls and other linear schemes. The secrecy loss is also shown as the gap between 
dashed and solid lines. 
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Fig. 5. Per-user secrecy rate vs. p for K = 4 users: without power allocation (solid), with «ls and p opt (dashed) and with joint 
optimal (a, p) op t (circle). The rate of the optimal RCI precoder without secrecy requirements (square) and the secrecy capacity of the 
MISOME channel (diamond) are also plotted. 
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TABLE I 

Proposed algorithms for power allocation. 



Algorithm 1 


Algorithm 2 


Initialize iteration counter t = 
Initialize all = 1, 6^ = 
repeat 

Solve l[56j to obtain p (t) 

Update 6^ at z = SINR fc (p (t) ) 

Increment t 
until convergence 
Obtain p k = e ?fc , fc = 1, . . . , if 


Initialize iteration counter ti — 0, <2 = 

Initialize pfc = 1/7, and set p k = logpfc, fe = 1, . . . ,K 

Initialize Qo = if Copt using equation d30t 

repeat 

Increment ti 

Obtain a| using steepest descent with a tl -i as initial point 

Initialize all a<* 2) = 1, = 

repeat 

Solve $56[ to obtain p ( ' 2) 

Update 4* 2) , &fc 2) at = SINR fc (a^ , p ( ' 2 ') 
Increment ti 
until convergence 
Set p = p (t2) 
until convergence 
Obtain p k = e? k , k = l,...,K 



